Diffusion-annihilation processes in weighted 
scale-free networks with identical degree sequence 

Yichao Zhang 

Department of Computer Science and Technology, Tongji University, 4800 Cao'an 
Road, 201804, Shanghai, China 

Zhongzhi Zhang 

E-mail: zhangzzOf udan . edu . cn 

School of Computer Science, Fudan University, 200433, Shanghai, China 
Shanghai Key Lab of Intelligent Information Processing, Fudan University, 200433, 
Shanghai, China 

Jihong Guan 

E-mail: jhguan@tongji.edu.cn 

Department of Computer Science and Technology, Tongji University, 4800 Cao'an 
Road, 201804, Shanghai, China 

Shuigeng Zhou 

E-mail: sgzhouOf udan . edu . cn 

School of Computer Science, Fudan University, 200433, Shanghai, China 
Shanghai Key Lab of Intelligent Information Processing, Fudan University, 200433, 
Shanghai, China 

Abstract. The studies based on A + A — > and A + B — > diffusion-annihilation 
processes have so far been studied on weighted uncorrelated scale-free networks and 
fractal scale-free networks. In the previous reports, it is widely accepted that the 
segregation of particles in the processes is introduced by the fractal structure. In 
this paper, we study these processes on a family of weighted scale-free networks 
with identical degree sequence. We find that the depletion zone and segregation are 
essentially caused by the disassortative mixing, namely, high-degree nodes tend to 
connect with low-degree nodes. Their influence on the processes is governed by the 
correlation between the weight and degree. Our finding suggests both the weight and 
degree distribution don't suffice to characterize the diffusion-annihilation processes on 
weighted scale-free networks. 



1. Introduction 



Complex networks are a powerful and versatile mathematical tool for representing 
and modeling the the structure of complex systems [HE]. Their wide applications in 



Diffusion- annihilation processes in weighted scale-free networks with identical degree sequenced 



distinct areas have made them an extensively focused discipline in the past decade [31 S]. 
Prompted by data mining and the increased computing power of computers, extensive 
empirical studies have unveiled that most real networked systems can be characterized 
by a power-law degree distribution P(k) ~ A; -7 , leading to the rising of research on 
our basic understanding of the organization of many real-world systems in nature and 
society jH [2j [31 H] . The characteristic exponent 7, usually observed in the range G (2, 3] 
in recent empirical studies (TJ [2j [31 0] is very important since it fundamentally influences 
some dynamical processes on the scale-free networks, e.g., synchronization 016], disease 
spreading [7], and so forth. Among these processes, one aspect that has recently received 
considerable attention is the diffusion-annihilation problem, i.e., bimolecular chemical 
and physical reactions of the identical particles A + A — > and different particles 

A + b ^ <b \®UM\n\ti%U3M\IR\I®- 

Unlike diffusion-reaction, the substances in diffusion-annihilation process don't 
yield products with mass. In the study of diffusion-annihilation, the density p of the 
surviving particles is thus a crucial problem since it presents a quantitative description 
of the reaction process. In the large time limit, p behaves as 
1 1 

^T-^T = k-t f , (1) 

P(t) P(0) 1 ; 

where k is the rate constant and p(0) is the particle density at t — 0. In the 
mean- field approximation with pa(0) = Pb(0), both processes can be described as 
= —const ■ pit) 2 } whose solution is / = 1. The solution is valid in regular 
lattices of Euclidean space [12] with a spatial dimension d > d c , where d c is the critical 
dimension of this process. For the A + A — > process, d c = 2 while for the A + B — > 
process d c = 4 [13]. Further studies on fractals found that the exponent / = for 
A + B — )• [14|, where d s is the spectral dimension of the fractal structure. 

As the existence of the depletion zone (A + A — > 0) [15] and segregation of the 
reactants (A + B — > 0) [16], the upper bound of the exponent / for the regular lattices 
is 1. Whereas, when the processes are performed on scale-free networks with identical 
nodes and links, / can be considerably higher than 1 [Sj. Inspired by the observations, 
the relation between 7 and / on A + A — > was investigated analytically [H] in 
uncorrelated scale- free networks [TTJ, [18] . Put briefly, the term "uncorrelated" denotes 
that no degree-degree correlations among nodes exist in the networks, namely, the 
conditional probability P(k'\k) that a node of degree k is connected to a node of degree 
k' can be formalized as k ^ ■* . The analytical solution shows / is only governed by the 
exponent 7 for this class of scale-free networks. Subsequently, an interesting study of 
A + B — > on fractal scale-free networks shows the segregation can also be found in the 
scale-free networks [10J. Influenced by the segregation, the reaction process is hampered 
apparently. 

Very recently, considering heterogeneous distributions of weights [191 120] , a heuristic 
research on the weighted uncorrelated scale-free networks analytically present a more 
realistic conclusion [11]. In this work, the weight of links is defined as Wij — (/cjfcj) 
with the degree hi and kj of both nodes, where 9 is the network's weightiness parameter 
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which characterizes the dependence between link weight and the node degrees [19]. 
When 8 = 0, there is no dependence between link weight and node degree, all link 
weights are equal with one, and the network becomes an unweighted network. When 
6 > 0, it is a weighted network where links have different weights. The larger 9 and the 
wider difference between links. Based on the mean-field rate equation for the average 
density pk of a node with degree k, for the A + A — > process, the authors showed 



in asymptotically large networks. For the A + B — > process, inserting the mapping 
relation, they claimed 



It has been shown that, / is only governed by the weight and degree distribution. 

In this paper, we study a family of weighted scale-free networks with the 
identical degree sequence (weighted IDS-SF networks), the reaction processes are vastly 
different from the previous reports [8j M, QUI [TT]. To this end, we briefly introduce 
weighted random diffusion in Section [2j Section [3] is devoted to explicit the IDS-SF 
networks. In Section HJ our extensive numerical simulations are compared with previous 
analytic results of the diffusion-annihilation processes running on top of the weighted 
uncorrelated scale- free networks [UJ . Finally, our conclusions are presented in Section [5j 
Our findings indicate that the disassortative mixing of the nodes is the essential reason 
for generation of the depletion zone and segregation in this class of scale-free networks. 

2. WEIGHTED RANDOM DIFFUSION 

Before introducing the construction of the networks, we briefly introduce the general 
random walk on weighted networks to clarify the influence of high-degree nodes (hubs) on 
the weighted random diffusion on scale-free networks. Random walk is a mathematical 
formalization of a trajectory that consists of taking successive random steps. A familiar 
example is the random walk phenomenon in a liquid or gas, known as Brownian 
motion [211 [22] . Random walk is also a fundamental dynamic process on complex 
networks [23]. Random walk in networks has many practical applications, such as 
navigation and search of information on the World Wide Web and routing on the 

Internet [2H E3 E51 EH EB] • 

Let's consider a weighted random walker starting from node % at step t = and 
denote Pi m (t) as the probability of finding the walker at node m at step t. The 
probability of finding the walker at node j at the next step is 




(2) 




(3) 




(4) 



rn 
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where element of the network's adjacent matrix. 

In this case, we define the weight of a link between nodes i and j as 

{0 link i-j doesn't exist 
(kikj) e link i-j exists 

where fc, and kj denote the degree of node i and j respectively. On the other hand, the 
strength of node i is defined as 

Si = ^2 w » = ( k i k i) 6 ' ( 6 ) 
j€T(i) jer(») 

Thus the probability Pij(t) for the walker to travel from node % to node j in t steps is 



x 



mi,...,mt-i 
W 



x ... X r^zlL, (7) 

s mt-i 

In other words, P^t) = J2 mi ,...,m t -i P i™x P m x mi • • • Pm t -ij- Comparing the expressions 
for Py and Pji one can see that SjPy(t) = SjPji(t). This is a direct consequence of the 
undirectedness of the network. For the stationary solution, one obtains P°° = Si/Z with 
Z = £\ Si. Note the stationary distribution is, up to normalization, equal to Sj, the 
strength of the node %. This means the higher strength a node has, the more frequently 
it tends to be visited by a walker. Notably, for degree uncorrelated networks [29], Sj in 
the steady state scales with ki as s« ~ kf +l [20] . 

3. The scale-free networks with identical degree sequence (IDS-SF 
networks) 

The scale-free networks with identical degree sequence are a common topic in complex 
networks, which offer researchers a platform to understand how the dynamical behaviors 
are influenced by the degree heterogeneity of networks [30j [31] . As a class of the these 
networks [321 [33] , the construction of the present model is controlled by a parameter 
q [32] [33] as shown in Fig. [H evolving in a recursive way. We denote the network after 
n iterations by G(n), n > 0. Then the networks are constructed as follows. For n = 0, 
the initial network G(0) consists of two nodes connected to each other by a link. For 
t > 1, G{n) is obtained from G(n— 1). That is to say, to obtain G(n), one can add three 
links to each link existing in G(n — 1) (as shown on the left of Fig. [T]) with probability 
q, or replace it with a quadrangle (as shown on the right of Fig. [TJ) with complementary 
probability 1 — q. In Fig. [2} next, we present the first three iterations of two special 
networks corresponding to two limiting cases q = and q = 1, respectively. 

As discussed in the reference [32], these two limiting cases and the middle cases 
(0 < q < 1) exhibit many interesting properties. For instance, the same degree sequence 
independent of parameter q, the identical degree distributions, and no triangles [31] 
formed by connections among the neighbors. Note that, as shown in Fig. [3] the Pearson 
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Figure 1. (Color online) Iterative method of the network construction. Each edge is 
replaced by the connected clusters on the left-hand side with a certain probability q, 
otherwise by the one on the right-hand side, where red squares represent new nodes. 
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Figure 2. (Color online) Illustration of the first three iterations of the network for 
the particular cases q = and q = 1. 



coefficient increases with q generally, indicating the IDS-SF networks are disassortative 
for q = (the index tends to —0.5 as N — >■ oo [35J) and uncorrelated for q = 1 [36J. 
Hence, for q = 1, the topological structure of network satisfies the conditions of applying 
P(k\k') = kP(k)/(k) [TT] and mean field approximation well, which will be discussed 
in Section 14.11 in detail. Adopting several q values from to 1, one can generate 
various networks, for example, fractal (q = 0) and non-fractal (q = 1) networks. These 
particular features have the kinetics taking place upon the model be distinct from the 
well known results for other networks, e.g., the Barabasi- Albert (BA) graph [9j[37] and 
uncorrelated configuration networks |9j [31] . In the following, we will show a number of 
interesting behaviors of Diffusion-annihilation processes on the networks. 

4. Diffusion-annihilation processes on the weighted IDS-SF networks 

According to the conclusion on the weighted random diffusion operating in the weighted 
uncorrelated scale-free networks in Section EJ it is easily seen that P°° = ' , e+1 . Thus, 

1 i 

for 6 > 0, particles move towards hubs with time gradually. As hubs are the minority 
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Figure 3. Pearson correlation coefficient, average path length, and diameter versus 
q ranging from to 1 for the IDS-SF networks. Each data point corresponds to ten 
independent realizations of the network for n = 7. 



of the population, moving to them means getting concentrated actually. At these hubs, 
particles have thus a high probability to collide and react with each other, leading to a 
higher reaction rate than that in homogeneous networks [9]. For 9 < 0, conversely, the 
particles are repelled by the hubs. In this case, the particles are getting dispersed on the 
low-degree nodes with time, which are also called leaves. Seen in this light, reaction rate 
of diffusion-annihilation tends to decrease with 9 for the two processes. In what follows, 
we will show the diffusion tendency mentioned above is correct, but the influence of 9 
on the reaction rate is not monotonic, which depends not only on degree distribution 
but also other topological features of the network. 

We first generate a special IDS-SF structure through an iterative way with n = 7. 
The simulation results are obtained on IDS-SF networks with 10, 924 nodes and 16, 384 
links. For the two reaction processes, each node in the networks can host at most one 
particle. The concrete processes are defined as follows: an arbitrary particle jumps with 
a certain probability — from a node i to a randomly chosen nearest neighbor j. If it is 
empty, the particle fills it, leaving i empty. If j is occupied, the two particles annihilate, 
leaving both nodes empty. An initial fraction p(0) of nodes in the networks is randomly 
chosen, which is occupied by an A particle with probability 0.5 for both types. For 
the A + B — > process, the initial densities of A and B are equal, i.e., ps(0) = Pa(0). 
For a convenience of discussion, we define / as the first order derivative of ^y, where 
p(t) = p A (t) for A + A -)> and p(t) = p A (t) + p B (t) for the A + B -)> process. 
In the cases among < q < 1, each plot corresponds to 100 simulations that are ten 
runs for ten independent realizations of the network with the same parameters. For the 
two limiting cases q = 0, 1, each plot corresponds to 100 runs for the two deterministic 
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networks. 

As the degree sequences of the IDS-SF networks with q G [0, 1] are the same, in 
which 7 = 3 [32], Eq. (J2J) and Eq. ([3]) can be rewritten as 

( t 1 e < o 

p ~ < o < e < i . (8) 

[ e~* > 1 

For the A + B — >• process, inserting 7 = 3, one can also obtain 

( r 1 e < o 

p~ J (tint)- 1 = . (9) 

[ O<0<1 

For a convenience, we define two quantities as follow: 

N AA (t) nm 

iVAfl(t) n , 

M(t)(M(t)-l)' 1 j 

where iV^^) denotes the number of close contacts between two nodes with the identical 
particles for the A + A — > process. NABif) denotes the number of contacts between 
the distinct particles for the A + B — > process at time t [38]. M(t) denotes the total 
number of particles at time t. 

4-1. Case of q = 1 

As shown in Fig. [2j in the case q = 1, the networks are reduced to the (l,3)-flower 
proposed in the reference [35]. By definition [39], the fractal dimension df can be 
obtained by 

d;=lim n^A (12) 



?1— >00 



In L 



where N n and l n are the size and diameter of G n respectively. Inserting N n = |(4 n + 2) 
and l n = 2n [35] into Eq. ( ITS]) , we have 

cL = lim (31n2n) . (13) 

n— >oo 

Obviously, the net is infinite-dimensional, namely, a non-fractal network. For A+A — > 0, 

Fig. H] shows the relation between / and time t for 9 = —1, 0, 1, where red lines are the 

power-law fittings of the plots. Note that all the plots about the dependence of / and 

i i 

Qaa(ab) on time t are logarithmically binned in this paper. Concretely, f t2 = p(i ^_ f P i (tl) 
and Qt 2 = ; where the time interval log{t2) — log(ti) = 0.1. For each panel, 

the curve is relatively stable in the beginning and fluctuates radically in the end. The 
numerical results show the reaction processes are vastly different from the previous 
analytical predictions on the weighted uncorrelated networks denoted by the dashed 
lines [TTj . 
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Figure 4. (Color online) / and Qaa as a function of time t for A + A — > with 
q = 1. The first and second row denote / and Qaa versus t respectively. The dashed 
lines correspond to the mean field prediction: For the A + A — > process, / = 1 when 
9 = 0, — 1, and / = exp(t) when 9 = 1. 



Compared with the mean-field prediction, one can observe many discrepancies in 
Fig. HI Here, we only focus on looking for some common reasons. Note that, for 9 > 1, 
the mean-field prediction / = e* is much higher than our results in Fig.H^a). To show the 
plots clearly, we omit the dashed line in this panel. The apparent discrepancy is mainly 
caused by the approximation N g — > oo. In this condition, (k 1+e ) — > oo, which makes the 
differential equation solvable. For < 9 < 1, on the other hand, the approximation in 
the literature [11] omits the reaction running on the low-degree nodes, which causes the 
predicted reaction rate is lower than our observation. For 9 = 0, our results in Fig. H(b) 
roughly match the conclusion in the reference [UJ in term of the scaling of /. But, the 
value of / is a bit higher than the prediction in that the global mean first-passage time 
of random walks G in the mean field prediction of Eq. (JS} is proportional to N g [ID] 

ln3 

while G ~ N g ln4 [33]. So that, one can expect a larger deviation in the IDS-SF networks 
with n > 7. Notably, this observation in this case is inconsistent with the previous 
conclusion on finite size effects, i.e., ~ N^~t for 7 < 3 [9J. For 9 < 0, our results 
in Fig. 11(c) are basically lower than the prediction. This deviation is caused by the 
approximation in Taylor expansion. As is known, p(t) can only be omitted at the end 
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Figure 5. (Color online) / and Qab as a function of time t for A + B — > with q = 1. 
The first and second row denote / and Qab versus t respectively. The dashed lines 
correspond to the mean field prediction: For A + B— >0, f = I when 9 = — 1 and 
/ = ln(t) + 1 when 9 = 0. 

of the reaction, where it is close to 0. 

For the A + B — > process, we measure the relation between the total particle 
density p(t) = pa{1) + Ab(^) and time t as shown in Fig. El Our observation exhibits the 
similar behaviors with A + A — > 0. As the probability of collision between two identical 
particles is equal to that for distinct ones, one can find Qab in Fig. El is about half of 
the corresponding Qaa- Thus, the reaction rate of the A + B — > process is naturally 
much lower than that of A + A — > 0. 

It should be mentioned that the (l,3)-flower is a network with a number of 
common properties, e.g., non-fractal topology, no degree correlations and scale-free 
degree distribution, satisfying the condition of mean- field approximation fully. But, 
the annihilation dynamics on it presents many unpredictable properties. Thus, there is 
a need to provide such a complement to the previous discussion on both fractal scale-free 
networks and weighted scale- free networks [TT]. Without loss of generality, in what 
follows, we will investigate the other limiting case q = 0. 
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Figure 6. (Color online) / and Qaa as a function of time t for A + A — > with q = 0. 
The first and second row denote / and Qaa versus t respectively. 

4.2. Case ofq = 

Unlike the case of q — 1 addressed above, for q = 0, the networks are reduced to the 
(2, 2)-flower as shown in the corresponding panel of Fig. [2J which is a fractal network 
whose fractal dimension dj = = 2 [35]. By definition, the fractal network is a network 
satisfying the fractal scaling Nb{Ib) ~ lj > "where Nb is the number of boxes needed to 
cover the entire network with boxes of size Ib- Note that the fractal scaling df holds 
in the system where hubs are located separately from each other [H] 02]. As is known, 
the mean-field theory can only be applicable when the nets have infinite dimensionality 
but not in the fractal ones [35}, \TU\. Thus, the discrepancies between the mean field 
prediction and our results are not unexpected. However, the weighted networks have 
their unique subtle properties, which gives rise to many interesting dynamical behaviors 
distinct from the previous unweighted fractal nets. 

For the A + A — > process in Fig. [6j we also measure the relation between / and 
t for the set of 9. In Fig. [61 one can observe that / decays with time t in all the three 
panels. In the case of 9 = 1 in Fig. [6](a), the exponent / decreases with t abnormally 
and exhibits a contrary behavior with the case of q = 1, in which / increases with t. 
For the A + B — > process shown in Fig. [TJ one can observe a similar phenomenon with 
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Figure 7. (Color online) / and Qab as a function of time t for A + B — > with q = 0. 
The hrst and second row denote / and Qab versus t respectively. 

the A + A — » process as well. Because of Qaa ~ 2Qab, f in this case is much lower 
than the A + A — > as well. Notably, for the unweighted case, i.e., 6 = 0, / shown in 
Fig. [3(b) is not a constant 0.5 mentioned in the reference |10j . 

For homogeneous initial distributions with equal densities of A and B, pa(0) = 
Pb(0); local hubs and the random fluctuation in the initial particle number generate the 
segregation of distinct particles, which drastically slows down the reaction rate. Usually, 
for the unweighted uncorrelated scale-free networks, it is hard for a large number of 
particles to form a close formation that cannot be penetrated by the other species 
because of a short diameter. However, for q = 0, the influence of disassortative mixing 
is enhanced by the high heterogeneous weight distribution as shown in Fig. [6(a) and 
Fig. [7(a). The tighter local hubs attract the particles, the lower the diffusion rate is. 

As shown in Fig. [8(a), a hub leads to a fast decay of the local A particle density in 
the beginning, followed by a slow decay in the long time regime as shown in Fig. [6(a). 
Thus, one can clearly observe depletion zones emerging from the intervals among hubs 
in this panel. In Fig. [6(b), a hub in a A or B-rich domain can give rise to a pure A or 
B zone after a prompt local annihilation of A and B, leaving a relatively particle-free 
space among the hubs. With these segregations, one can observe a slow decay of the 
reaction rate as shown Fig. [7(a). 
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Figure 8. (Color online) Illustration of the A + A — > and A + B — >• processes on 
the IDS-SF networks with n = 4. (a) A + A -> at t = 10, and (b) A + B -> at 
t = 15. Red and blue plots denote A and B particles respectively. 

Interestingly, the depletion zone and segregation also inhibit particles moving from 
leaves to hubs when 9 < 0. The behavior can be observed by measuring the average 
degree of occupied nodes increases with t as plotted in Fig. EJ Note that the plots are 
also logarithmically binned. Recalling the discussion at the beginning of this section, 
particles are attracted by the leaves in this condition. For q = 1, particles tend to 
agglomerate around hubs for 9 > (see Fig. [9(a)) and leaves for 9 < (see Fig. Etc)). 
Comparing Fig. [4(a) and Fig. [5(a) with Fig. [6(a) and Fig. 17(a), one can figure out the 
depletion zone and segregation slow down the rate of reaction for 9 > 0. For 9 < 0, 
comparing Fig. [5(c) with Fig. [7(c), instead, they accelerate the rate slightly. This is 
because they slow down the decentralizing process of the particles. For 9 = 0, their 
effect is hardly identified as shown in Fig. [9(b) without the enhancement of weight. 
Apparently, these interesting behaviors are vastly distinct from the previously reported 
results PUE]. 

4.3. Case of0<q<l 

For < q < 1, the networks are stochastic, which makes them not self-similar [33]. Thus 
the networks are non-fractal in this middle case. Thus, in order to discuss the variation 
in the dependence of / on q, we have performed extensive numerical simulations for 
various q from to 1. The simulation settings were the same as the former cases. When 
q increasing from to 1, the exponent of global mean first-passage time of random walks 
G(N g ), decreases from 1 to p| [33], which indicates the enhancement of transporting 
efficiency during the process. At the same time, the diameter of the networks also 
decreases while the disassortative mixing feature disappears. 

In Fig. CGU one can observe that the segregations among hubs disappear gradually 
with the increase of q. Under the influence, the diffusion rate increases drastically, 
leading to an apparent enhancement of / for 9 = 1 (see Fig. [10( a)). Notably, in panel 
(b), the purely topological segregations for 9 = are also observable, although its 
influence is not as apparent as that in panel (a). Also, comparing q = 0.1 with q = 0.9, 
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Figure 9. (Color online) (k) as a function of time t for A + A — > and A + B — > 
for = —1,0,1 with q = 0, 1 respectively. 



the subtle influence of segregations on the reaction rate in the case of 6 = — 1 can be 
identified in panel (c). As shown in this panel, the reaction rate decreases slightly with 
q, which is consistent with our observation in Section 14.21 . 

5. CONCLUSION 

In summary, we have investigated the diffusion-annihilation process on a family of 
weighted scale-free networks with identical degree sequence (weighted IDS-SF networks), 
which is controlled by a parameter q £ [0,1]. In this paper, the weight of links is 
defined as Wij = (hkj) 9 with the degree ki and kj of both nodes, where 6 is the 
network's weightiness parameter. For a convenience, we define a kinetic exponent / 

as —tfp-, where p(t) = pAif) for the A + A — > process and p(t) = pA(t) + 
for the A + B — > process. Based on the definition, we provide numerical results to 
characterize the relation between / and the reaction time t for the A + A — > and 
A + B — y bimolecular reactions. 

One significant observation is that, in contrast to the commonly accepted conception 
that the depletion zone and segregation only exist in fractal networks. Our observation 
shows they can exist in the diffusion-annihilation process on the non-fractal networks as 
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Figure 10. (Color online) Reaction rate exponent / with 6 = —1,0,1 and q = 
0.1,0.5,0.9 for A + A -> (black lines) and A + B — > (red lines) processes on the 
IDS-SF networks respectively. 

well. This striking feature in scale-free networks was not reported in previous studies. 
In fact, the depletion zone and segregation can both exists in fractal and non- fractal 
networks no matter whether it is weighted or not. We found that the segregation effect 
is essentially caused by the disassortative mixing, i.e., high-degree nodes tend to connect 
with low-degree nodes. On the weighted networks, its influence on the particles diffusion 
is highly enhanced by the weight heterogeneity. We have demonstrated that both degree 
and weight distribution do not suffice to characterize the diffusion-annihilation processes 
on weighted scale-free networks. Our observations suggest care should be taken when 
making general statements about the diffusion-annihilation process in weighted scale- 
free networks. 
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